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ONLINE WAVEERONTMEASU^ 

BACKGROUND OF THE INFORMATION 

Field of the Invention . 

Embodiments of the invention are generally directed to the field of 
image processing, and more particularly to methods and apparatus for the 
simultaneous measurement, analysis, and display of ocular wavefront 
infomiation referred to herein as "online*' aberrometry. 
Description of Related Art 

There are few technologies that have not been touched by the science 
of image processing. Advancing developnient in mpst fields typically gives 
rise to process automation where the ability to quickly and accurately detect 
structures in camera images becomes increasingly important One area 
targeted by the embodiments of the instant invention is wavefront aberration 
measurements of the human eye, although the techniques set forth herdn 
below will find much wider application directed to the detection of structure 
in images. 

The field of ophthalmology concerned with vision correction through 
refiractive surgery, or the iise of lens components on or in the eye, has 
experienced the relatively recent development of aberrometers. These 
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devices allow practitioners to explore the meaning of vision quality and to 
appreciatethatitisnotlimitedmerely to visual acuity. Aberrometers which 
incorporate wavefront sensors can provide information about vision defects 
that, upon correction, may not only provide visual acuity at the theoretical 
limit but also better vision, perhaps even customized vision, under a panoply 
of viewing conditions. 

One of flie most conventional and well developed forms of ocular 
wavefront sensing relies on the Hartmann-Shack principle. A 
Hartman-Shack wavefront sensor typically includes a microlens array that 
iinages various portions of a distorted wavefront exiting the eye onto a CCD 
detector/camera. The image produced by the microlens array con^rises an 
array of small dots of light that are sUghtly displaced from reference 
locations of the light dot image from an unaberrated wavefront The 
aberrated dot displacements are related to the localized slopes of the 
wavefront exiting the eye's pupil. Zemike polynomials (or other 
mathematical forms) can be derived from these displacements, which are 
then used to characterize virtually all of the eye's aberrations. The ability to 
make accurate wavefront calculations is critically dependent upon the trae. 
determination of the center location of each dot in the wavefront ima ge. 
This aspect of the wavefront analysis process is known as centroid 
detection. 
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Hartmann-Shack wavejfront sensors, and otiier well-known types such 
as Tscheming, for example, typically measure single images of centroids or, 
at best, a very small number ofimages over a short time interval. The eye, 
however, is a dynamic system with rapidly varying wavefront changes. The 
time needed for centroid detection has been the primary culprit hindering 
real-time measurements with repetition rates greater than a few.iniages per 
second. A system kaown in the wavefront art as WASCA has dernonstrated 
a repetition rate of about 7Hz for the 30-sec6iid record of a wavefront 
However, the wavefront images must first be recorded, saved, and 
subsequently evaluated. A single wavefront image requires about 400Kb of 
computer memory. Moreover, aberration measurements (e.g,, sphere, 
cylinder/axis, and higher order aberrations) cannot be displayed online, i.e., 
substantially simultaneously with the wavefront measurement and 
calculation. Nor is it possible to acquire and save pupil images and centroid 
images substantially simultaneously, making it virtually inrpbssible to 
evaluate the influence of eye movement on changes in the wavefront These 
illustrations represent some of the exemplary development issues in the field 
of ocular wavefront measurement addressed by the embodiments of the 
instant invention. 
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Brief bescription of the Drawings 

Figure 1 is a flow chart schematic of a fast centroid detection 
algorithm embodiment of the invention; 

Figure 2 is a reproduced photographic of a Hartmann-Shack 
wavefront (centroid) image associated with an embodiment of the invention; 

Figure 3 is a drawing representation of a wavefront (centroid) image 
associated with an embodiment of the invention; 

Figure 4 is a flow chart schematic of a sorting algorithm according to 
an embodiment of the invention; 

Figure 5 is another drawing representation of a wavefront (centroid) 
image representing an algorithm process according to an embodiment of the 
invention; 

Figure 6 is a schematic illustration representing another algorithm 
process according to an embodiment of the invention; 

Figure 7 is an illustration representing another algorithm process 
according to an embodiment of the invention; 

Figure 8 is a reproduced photographic of part of a Hartmann-Shack 
wavefront (centroid) image representing another algorithm process 
according to an embodiment of the invention; 

Figure 9 is a flow chart schematic directed to another embodiment of 
the iQvention; and 

4 
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Figure 1 0 is a block diagram of an apparatus embodiment of the 
invention. 

Detailed Description of Preferred Embodiments of the Invention 
An embodiment of the present invention is directed to an algorithm 
that provides for the detection of bright and dark stractures in a CCD image. 
The image may comprise, for illustration pmposes only, centroids in a 
wavefront image, or identified markings on a contact lens surfece. Figure 3 
illustrates an array of centroid images 32 as black spots in a simulated 
Hartmann-Shack wavefiront image. 

A feature of this embodiment is that the positions of stractures such 
as centroids 32 in a CCD image, for example,-can be located and sorted in 
approximately 5ms, Zemike coefficients calculated, and wavefront image 
information displayed in approximately 13ms or less (utilizing an 8OOMH2 
Pentium* (Intel) processor or equivalent). A 1.6GEz processor would 
reduce the total time by approximately one-half by increasing the rate to 
approximately 50Hz. The algorithm thus describes what is referred to 
throughout this description as "online" analysis of wavefront information, 
referring to the detection, measurement, and display of wavefront 
. information substantially simultaneously at a rate up to approximately 
25Hz. • 

5 
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The embodimeiits to be described in detail below will be illustrated in 
terms of a Hartmaxm-Shack wavefront measurement and analysis, but .are 
more broadly applicable to other wavefront methodologies, for example, 
those based on the Tscheming principle and others known in the art; 
however, as noted above, the embodied invention is even more generally 
applicable to image processing for analyzing stmcture in a CCD type image. 

A method for centroid detection in a wayefront image according to 
the instant embodiment is set forth in the flow chart diagram 100 of 
Figure 1. It is noted at the outset that in the illustrated embodiment, a 
sequential plurality of images are acquired, analyzed, and displayed as 
desired at a rate of 25hz, but for the sake of simplicity the algorithm steps 
set forth below apply to a single wavefront image and are repeated for every 
wavefront image of the desired sequence. 

In Figure 1 at step 102, an XxYpixel size wavefront image, as 

shown for example by image 22 in Figure 2, is acquired. The light spot 

images are represented by variable pixel signal intensities as shown. Images 

taken from a CCD camera usually consist of an array of pixels in which 

every pixel is assigned a number proportional to the amount of charge 

collected at this pixel. This number is referred to as the signal of the pixel. 

In the illustrative description that follows, a regular square grid of bright 

points in a dark image will be described in detail. 

6 
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i) Image Compression 

After acquiring the image at step 1 02, the image is compressed from 
size X X Y pixels to X/n x Y/m pixels at step 104. This can be done by 
averaging the signal for every pixel in an n x m pixel square in the original 
image starting, for example, m the upper left comer of the image and 
scanning through the iniage. The signal in the upper left comer of the 
compressed image is then set to the average pf the first square, the signal of 
the next pixel is set to the average of the next (second) square, and so on, 
finally yielding a picture X/n x Y/m pixels in size., n and m should be 
integers with X/n and Y/m also being integer values. In an.exen]plaiy 
embodiment n = m = 8. 

ii) Background Subtraction 

At step 1 06, the compressed image is then divided into square 
regions, or tiles (not to be conftised with the pixel squares in (i) above). In 
the exemplary embodiment, one tile is a square of 64x64 pixels, but other 
sizes can be used. Typically, one tile might contain 3-5 centroids. The 
average signal is again calculated for every tile. The average values for the 
tiles are then linearly extrapolated to yield a background value for any 
location in the image. This background is then subtracted from the image 
yielding a low signal outside the centroids. In the illustrative embodiment, a 
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signal-to-noise ratio of two was improved to a signal-to-noise ratio of 10 by 
the background subtraction. 

iii) Rough Structure Detection 

At step 108, the approximate, or rou^ stracture points (centcoids) 
are identified. First, a maximum is defined as the highest signal in the 
compressed image. The maximum is determined by a scan through the 
image, and the X-position, Y-position, and signal value of every pixel is 
recorded hi a table, but only if the signal value of this point is greater than a 
certain percentage of the maximum, e.g., 30 % (other values can be selected 
by the user). In the exemplary embodiment, this yields a table (Table I) of 
about 400 entries. The table is sorted by descending signal as shown. Any 
of a variety of quick sorting routines is available to accomplish this. 

TABLE I 



Signal 


X Position 


Y Position 


223 


86 


55 


154 


85 


75 


135 


87 


95 


133 


115 


56 


110 


118 


74 


108 


114 


93 


• 


■ • 


• 


• 


• 


• 


« 


• 


• 
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The first entry (highest signal value) is defined as the first rough structure 
point. Then, all entries in the table that obey a certain pre-set condition are 
defined as rough structure points. In the exemplary embodiment, the pre-set 
condition is that the position of the particular entry is farther away fi-om all 
yet found rough structure points than a pre-set distance. In an exemplary 
embodiment, the distance is 17 pixels. After this first iteration, a table of 
rougih structure points is created that includes approximately 95% of all 
points that are to be detected. 

iv) Refine Detection of Structure 

To increases the confidence level that all points of the structure are 
found, step 108 can be repeated as shown at block 110, setting a new 
.min imu m to a certain percentage of the minimum in the first iteration. The 
second iteration finds points that were too weak in signal to be found in the 
first iteration. The rough structure points found in the first iteration are • 
accoimted for so they will not be found again (i.e., they do not obey the 
condition of being farther away than a pre-set distance fi:om the detected 
points). 

v) Ultimate Structure Detection • 

At step 1 12, the ultimate centroid positions are deteimined. Since, the 
image was earlier compressed in step 104, much of the information ' 

originally contained in the image was ignored. This information can now be 

9 ■ 
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used to determine more exact centroid positions. Using the original 
uncompressed image, a square of, for example, 15x15 pixels is created 
arbimd every rough point. Generally, each square is smaller than 2x the 
Tnim'mum distance to insure that each square contains only one centroi^ and. 
is larger than the centroid itself, la the exemplary embodiment this value is 
between five and 3 0 pixels. Then, the center of mass of signal for the signal 
distribution inside the square is determined, giving rise to the substantially 
exact position of the centroid. 

In an aspect of the embodiment, step 112 can be repeated, for 
example, 1 , 2, 3 . . .n times to determine still more exact results. The 
calculated center of mass in the previous step is subsequently used. Each • 
stmcture point can also be assigned a quality factor depending on how much • 
the position of the center of mass changes if the square around the pixel is 
willingly shifted by a user-set distance, at step 112. In an esemplaiy 
embodiment, this distance is five pixels. The points whose positions have 
changed the least are assigned the highest quality factors. In this manner, 
spurious points or noise assigned a low quality factor can be eliminated, as 
they likely represent false structure points. 

In the illustrative embodiment directed to Hartmann-Shack wavefiront 

sensing, it is desirable to be able to correlate the centroids with the 

corresponding image fonning lenslet of the microlois array.- Thus, an 

10 
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aspect 200 of the embodiment as illustrated in Figure 4 is directed to the 
process of sorting the detected centroids so sis to assign them to a regular 
square grid pattern. It will be appreciated by those skilled in the art that the 
algorithm can be easily adapted to other structures or conJBgurations such as, 
. for example, points on rings, or any straight line of points. 

At step 202, the desired sorting configuration is selected. . In the 
exemplary embodiment, the configuration is a square grid based upon the 
geometry of the microlens array. For every previously found centroid point, 
i, the formula for a straight line i$ calculated containing the centroid point, i, 
and having a slope of 1 (45°), as shown at step 204. For starting positions 
of the upper left comer or lower right comer of the image, slope values 
between 6.1 to 0.9 can be used. Likewise, when the starting position is the 
upper right comer or the lower left comer of the image, slope values firom 
-0.1 to -0.9 can be selected. At st&p 206, the distance 502 (nO, between the 
line 514 and, in the illustrative embodiment, the upper left comer 506 of the 
image 510 is calculated, as illustrated in Fig. 5. All centroids, i, are then 
sorted by Ui at step 208 starting with the centroid having the smallest Ui 
value. At step 210, the centroid with the smallest nj value is assigned to 
Row 1 and is stored in memory as the last centroid of Row 1. In an aspect 
of the embodiment, the last centroids of the existing rows are stored in 

memory during step 210. At step 212, a region 610 is defined that, in an 

11 
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exemplary embodiment, comprises an area to the right of the last centroid 

605 of a given row, having dimensions that can be controlled and varied by 

parameters of the lenslet array, and having a shape that is suited to detecting 

the selected grid configuration, as illustrated in Figure 6, which shows the 

search area 610 for the next centroid. Any shape suited fot detecting other 

grid configurations is alternatively possible. Examples oif the lenslet iarray 

parameters include maximum angle 702, Tniniminn distance 704» mayinmiTn 

distance 1 (706), and maximum distance 2 (708), as illustrated in Figure 7. 

Then, at step 214, the next higher ni value is selected and that centroid is 

checked, with respect to all existing rows, whether the centroid is in the 

defined regioiL If yes, then at step 216, that centroid is assigned as the last 

c^troid of that row. If no, than that centroid is assigned the last centroid of 

a new row. Steps 214-216 are now repeated for all centroids. In this 

manner, rows start to build up from left to rigjht At step 220, the average y- 

position for each row is calculated and the rows are sorted according to their 

average y-position. This step facilitates marking of the top row as Row 1, 

the next row as Row 2, and so on. 

Before describing the stqps for sorting the columns, it is beneficial to 

point out that the situation can occur as illustrated by the faintly seen points 

lying along lines 802, 804 as shown in Figure 8; i.e., some centroids 812, 

8 14 in the middle of a row have not been detected due to bad quality of the 

12 
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centroid points, and the cenfroids to the left and right have been assigned to 
different rows. • In this event, optional step 222 involves merging th^e i»ws. 
This is accomplished by the following sub-steps: From the mean average y- 
position for each row jfrom step 214, calculate the mean distance between 
the rows by subtracting yrowi-yiow2 (yielding the distance between Rowl and 
Row2); yrow2-yiow3 (yielding the distance between Row2 and Row3); and so 
on, and then taking the mean values of the obtained distances. In tlie 
exemplary embodiment, the criteria for merging rows j and k is: If yj-yk < 
f*aand(Pk^>Pj.,a8t orPyart<Pitot), 
where: 

f is a variable parameter in the range between about 0. 1-0.7, that is set by 
the user. In the exemplary embodiment, vaiues between 0.3 to 0.5 are used; 
a is the mean distance between rows (see above); 

is the X value of the first (left-most) centroid of tbek row; and 
Pkjast is the X value of the last (right-most) CMitroid of the k row. 
In other words, the rows are merged if they are much closer in the y-position 
than typical and if they don't overlap i.e. row j is either completely to the 
left or completely to the right of row k. 

The process for sorting the columns begins at step 224 where the list 

of sorted centroids by distance value from step 208 is used again. The 

centroid with the smallest ni is assigned to Column 1 and is stored in 

13 
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memory as the last centroid of Column 1. In an exemplary aspect, the last 

centroids of an existing column are always stored in memory during step 

224. At step 226, a region is defined that, in the exemplary embodiment, 

comprises an area below the last centroid of a given column having 

dimensions and shape that are controlled and varied by the same parameters 

of the lenslet array as set forth above. This is illustrated by tilting the 

diagram in Figure 6 downward by 90 degrees. At step 228, the next higher 

Ui value is selected and that centroid is checked, with respect to all existing . 

colunms, whether the centroid is in the defined region. If yes, then at step 

230, that centroid is assigned as the last centroid of that column. If no, than 

that centroid is assigned the last centroid of a new column. Steps 228-230 

are now repeated for aU centroids. In this manner, columns start to build up 

fi:om top to bottom. At step 234, the average x-position for each column is 

calculated and the colunms are sorted according to their average x-position. 

This step facilitates marking of the left-most column as Column 1, the next 

column as Column 2, and so on. 

The situation can occur, as mentioned above with reference to Figure 

8, where some centroids in the middle of a column have not been detected 

due to bad quality of the centroid points, and thus the centroids above and 

below have been assigned to different columns. In this event, optional step 

236 involves merging these colunms. This is acconaplished by the 

14 
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following sub-steps: From the mean average x-position for each column 
from step 228, calculate the mean distance between the columns by 
subtracting XcohmmrXcoiunna (yielding the distance between Columnl and 
Column2); Xcoiunin2-Xcoiumii3 (yielding the distance between Column2 and 
Columns); and so on, and then taking the mean values of the obtained 
distances. In the exemplary embodiment, the criteria for merging colvrams j 
and k is: If Xj-Xk < fa and (Pk^ > Pj4a8t or Pk,i,gt < Pj^ ), 
where: 

f is a variable parameter in the range between about 0. 1-OJ, that is set by 
the user. In the exemplary embodiment, values between 0.3 to 0.5 are used 
a is the mean distance between columns; 

Picfiist is the y value of the first (top-most) centroid of the k column; and 

Pji;ia8t is the y value of the last (bottom-most) caitroid of the k column. 

In other words, the columns are merged if they are much closer in the x- 

position than typical and if they do not overlap; i.e. column j is either 

completely above or completely below column k. 

From the sorted centroid positions, a Zemike calculation can be made 

to determine the wavefront aberration. According to the illustrated 

embodim^t, the processes for acquiring a single image and displaying the 

wavefront aberration information, using an 800MHz Feathim processor, are 

Hsted below with their corresponding processing times: 

15 
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Finding and sorting the centroids: ~5ms; 

Perfonning the Zernike calculations: ~5ms; 

Imaging operations: ~8ms; 

Image display: ~8ms 

Pupil coordinate location (optional): ~ 6-8ins; 

Contact lens position (optional): ~5ms 
In an illustrative aspect of the embodimait, the Zernike calculation proems 
is performed twice during a measurement cycle, once for the second-order 
terms and once for the higher-order terms. The total time per image is 
approximately 40nis or slightly less at a repetition rate of 25Hz. In contrast, 
ccaiventional real-time (as opposed to online) wavefront anal^is consists of 
storing a sequence of images and subsequently analyzing the images for 
wavefront information. These techniques are limited by computer storage 
requirements. For example, storing two images at 25Hz jR)r30 seconds 
requires approximately 400MB/measurement On the other hand, according 
to the illustrative embodiment, tiie storage of images is not necessary 
because the information contained in those images has already been 
extracted by detecting the centroid and pupil position. Storing only pi^il 
and centroid position data results in data storage requirements of only about 
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. 700Kb of memoiy for a 20 second measurement at 25Hz, which will yield 



500 images. 



Another embodiment of the invention is directed to a device-readable 



the algorithm described above. Appropriate, media are well known and 
include^ without limitation, CD, DVD, diskette, hard drivej carrier wave, 
and others. 

Another embodiment of the invention is. directed to a method for 
wavefront analysis. The method is set forth in the block diagram of Figure 
9. The method comprises at step 902 acquiring a plurality of wavefront 
images of light exiting a pupil of an eye, where each of the images includes 
a displaced centroid that is indicative of wavefront measurement, 
information of the ^e, and at step 904 calculating and displaying tiie 
wavefront measurement information online for a selected aberration order. 
As used herdnabove and throughout the description, the term "online** 
refers to the substantially simultaneous measurement, analysis, and display 
of the wavefront measurement information. The exemplaiy fast algorithm 
described in the previous embodiment is but one way to facilitate an online 
process. Persons skilled in the art will appreciate that oilier algorithms can 
be developed, or are currently available, that will also facilitate the online 
technique. 



medium that has stored thereon an executable instruction for carrying out 
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In an aspect of this embodiment, the images are acqmred at a rate 
equal to or greater than lOHz. In an exemplary aspecl; the images are 
obtained at a rate of 25Hz; however, computer processor speed is the 
limiting factor in the image acquisition rate. Jn another aspect, at least 
50 sequential images are acquired. In another aspect, the step 906 of 
calculating and displaying the wavejfront measurement inforination online 
for a selected aberration order is performed for a selected pupil diameter. 
That is to say, if a pupil diameter value is not selected prior to the- . 
measurement, the display will be limited to second order aberrations (sphere 
or spherical equivalent, and cylinder/axis) because second order aberrations 
are pupil diameter independent However, one may want to insure that the 
resultant wavefront information is interpreted for a substantially constant 
pupil diameter, or a controlled variation in pupil diameter. On the other 
hand, a pupil diameter value can be selected prior to measurement allowing 
any Zemike order (e.g., coma, spherical aberration, higiher-orders) for that 
given diameter to be displayed online. An exemplary range of pupil 
diameters is between about 2rmn to 10mm. 

In an exemplary embodiment, at a repetition rate of 25Hz, 500 images 

can be obtained (and thus 500 njeasurements can be made) over a 20-second 

time interval. The first 125 images might be obtained for pupil diameter, D 

< the next 250 images for Dn±, < D < Dj^^ and the remaining 
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125 images for Dmax < D. Using the set of 250 images, Zemike amplitudes 
for Z2xx thorough Zioxx can be calcidated as set forth at ste^ If desired, at 
step 910 an average value of a selected Zemike order can be calculated and 
displayed online. If an average value has been calculated, well known 
statistical techniques can be used to determine a standard deviation at step 
912 that provides an error band for the average wavefrpnt measurement In 
an aspect of this embodiraent, blinking periods can be determined which 
contain anomalous wavefront information, and the information during these 
blinking periods discarded. Likewise, if contact lens position were being 
measured, for example, it would be advantageous, to eliminate 
measurements during blinking periods as lens settlement takes a short time 
to occur. Knowing the value of the most firequentily occurring wavefront 
amplitude for a particular aberration order based upon pupil diameter will 
allow the practitioner to prescribe a treatment or vision correction resulting 
in optimum vision for the patient 

In a related asjpect according to this embodiment, at step 914 a 
sequence of pupil images corresponding to the wavefront images can also be 
obtained. These images can be saved simultaneously so that the influence . 
ofeye movement on changes in the wavefront can be evaliiated . 

Another embodiment of the invention is directed to a wavefront 

measuring apparatus 1000 illustrated by the block diagram in Figure 10. 
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The apparatus includes an illumination component 1010 that delivers a 
small spot of light onto the retina of the eye 1020, an imaging component 
1030 that forms a centroid image of the illumination light scattered from the 
retina aad exiting the pupil of the eye, a detector 1040 adapted to acquire the 
centroid image, a processor 1050 working in cooperation with the detector 
that executes a centroid displacement calculation to determine the wavefront 
measurement information, a display component 1060 operatively connected 
to the processor that displays a selected wavefront measurement 
information, and means 1070 for instracting an online calculation and 
display of tiie selected wavefront measurement information. An optional 
pupil camera 1080 andpupilometer 1090 is also shown, where components 
1 0 1 5 are beain spUtters or optical equivalents thereof 
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We claim: 

1 . A method for wavefront analysis, comprising: 

a) acquiring a plurality of wavefront images of light exiting a 
pupil of an eye, wherein each of the images includbs a displaced 
centroid that is indicative of wavefront measurement information of 
the eye; and 

b) calculating and displaying the wavefront measurement 
information online for a selected aberration order. 

2. The method of claim 1 , comprising acquiring the images at a 
rate equal to or greater than lOhz. 

3. The method of claim 1 or 2, comprising acquiring at least 50 
sequential images. 

4. The method of any of claims 1 to 3, further comprising 
determining an average value and a corresponding standard deviation 
for the selected aberration order wavefront measurement information. 

5. The method of any of claims 1 to 4, wherein step (b) is 
performed for a selected pupil diameter. 

6. The method of any of claims 1 to 5, comprising displaying the 
average value for the wavefront measurement information for a 
selected pupil diameter. 
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7. The method of any of claims 1 to 6, wherein the wavefront 
measiirement information is a manifest refraction value. 

8. The method of claim 7, comprising displaying the manifest 
refraction value online. 

9. The method of any of claims 1 to 8, wherein the wavefront 
measurement information is selected to correspond to one of the 
second through the 1 0* Zemike aberration orders or optical 
equivalents thereof. 

10. The method of any of clauns 1 to 9, comprising acquiring a 
correspondmg pupil image of the eye for each wavefront image; 

determining at least one of position data, size data, shape data, 
and geometric characterization data of the pupil in each pupil image; 
and 

displaying at least one of the pupil images and the 
corresponding pupil image data online. 

1 1. The method of claim 10, comprising displaying a selected 
wavefront aberration online, selected from a group comprising the 
second through the 10* Zemike aberration orders or optical 
equivalents thereof. 

12. The method of claim 10 or 11, comprising saving a pupil image 

and a temporally corresponding wavefront image simultaneously. 
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13. The method of claim 5, wherein the selected pupil diameter is 
in a range between about 2mm to lOmnL 

14. The method of any of claims 1 to 13, comprising calculating 
the wavefront measurement information utilizing the algorithm of any 
of claims 15 to 35. 

15. An algorithm for centroid detection in an image, comprising: 

a) acquiring an XxY size image represented by a variable pixel 
signal intensity; 

b) compressing the XxY size image to an X/n x Y/m size 
image, where n, m equal any integers and X/n, Y/m are integer 
values; 

c) determining a background uitensity for any position in the 
compressed iipage and subtracting this background from the 
compressed image; 

d) detecting a plurality of approximately positioned centroids in 
the background-subtracted compressed image; 

e) iterating step (d) until approximate positions of k desired 
plurality of centroids are detected; 

f) converting the approximate position of the desired plurality 

of centroids into more exact positions in the XxY size image, 

whereby every centroid position in the image has been identified. 
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16. The algorithm of claim 15, comprising: 

g) iterating step (f) until a desired level of more exact positions 
is determined. 

17. The algorithm of claim 16, furflier comprising: 
assigning a quality factor to each centroid in relation to a 

magnitude of positional change for each centroid in each iteration of 
step(g). 

18. The algorithm of any of claims 15 to 17, conq)rising: 
sorting the centroids determined from step (f) according to a 

predetermined configuration. 

19. The algorithm of claim 18, wherein the configuration is a 
geometric grid. 

20. The algorithm of claim 19, comprising a rectangular grid. 

21 . The algorithm of claim 1 8, wherein the configuration is a ring. 

22. The algorithm of claim 18, wherein the configuration is a 
straight line. 

23. The algorithm of claim 18, comprising: 

associating each determined centroid with a respective centroid 
image forming element 

24. The algorithm of any of claims 15 to 23, wherein compressing 

the XxY size image to an X/n x Y/m size image comprises: 
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a) averaging the signal for every pixel in an n x m square . 
starting in a furst predetermined region of the original image and 
scanning through the image, setting a signal level in a corresponding 
jfirst predetermined region of the compressed image to the average 
value of the first predetermined region; 

b) repeating step (a) for a second and subsequent 
predetermined regions until the X/n x Y/m image size is obtained. 

25. The algorithm of claim 24, wherein n = m = 8, 

26. The algorithm of claim 24 or 25, wherein the first 
predetermined region is the upper left comer of the image. 

27. The algorithm of any of claims 15 to 26, wherein step (c) 
con^rises: 

dividing the compressed image into a plurality of image 
segments each of which contains a plurality of centroids, determining 
an average signal value for each image segment, and extrapolating the 
average values for each image segment to determine the background 
intensity level. 

28. The algorithm of claim 27, wherein the image segriiehts are 
squares. 

29. The algorithm of claim 27, wherein each image segment 

contains approximately 3 to 5 centroids. 
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30. The algorithm of any of clauns 15 to 29, wherein step (d) 
comprises: 

a) detenmning a maximum signal value in the image; 

b) setting a threshold value as a predetermined percentage of 
the maximum; 

c) determining an X-position, a Y-position, and a signal 
strength of each pixel that has a signal strength greater than the 
threshold value; 

d) sorting the values from step (c) in descending order of signal 
strength; 

e) assigning the highest signal strength as first approximately 
positioned centroid; and 

f) selecting a pre-set condition for defining all sorted values as 
approximately positioned centroids, which obey the pre-set condition. 

31. The algorithm of claim 30, wherein the pre-set condition is that 
the position of each subsequent approximately positioned centroid is 

a farther distance away than a pre-set distance from all yet determined 
approximately positioned centroids. 

32. The algorithm of claim 30 or 3 1, further comprising setting a 

new threshold value to a predetermined percentage of a minimnm 

value of the sorted signal strengths and iterating steps (c-f), wherein 
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the already identified approximately positioned centroids are not 
identified again. 

33. The algorithm of any of claims 15 to 32, wherein step (f) 
comprises; 

defining a boundary structure around every approximate 
position of the desired plurality of centroids in the original image; and 

determining a center of mass of the signal for the signal 
distribution inside of the boundary. 

34. The algorithm of claim 20, con]{)rising: 

a) calculating a straight line formula for each sorted centroid, i, 
containing the centroid point, i, and having a slope between the values 
of about <-0.1 or >0.1; 

b) calculating a distance, Ui, between the line and a reference 
position in the image; 

c) sorting all centroids, in, by Ui starting widi the smallest oi 

value; 

d) assigning the centroid with the smallest Ui to a first row and 
storing this centroid as a last centroid in the first row; 

e) defining a reg;ion as ah area to the right of a last centroid of a 
given row having dimensions that are variably controllable by an 
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imaging component parameter and a shape suitable for detection of a 
selected grid structure; 

f) obtaining the next Ui value and determining, for all existing 
rows, whether the centroid is within the region; 

g) assigning the centroid as the last centroid in the given row if 
the centroid is witibin the region, or, assigning the centroid as the last 
centroid in a new row if the centroid is outside the regidn; 

h) repeatmg steps (f-g) for all centroids; 

i) calculating an average y-position for each row and sorting 
the rows according to the average y-positions to identify a top row. 
Row 1, a next row. Row 2, and so on to Row n; 

j) assigning the centroid with the smallest Uito a first column 

and storing this centroid as a last centroid in the first column; 

k) defijiing a region as an area below the last centroid of a 

given column having dimensions that are variably controllable by the 

imaging component parameter and a shape suitable for detection of 

the selected grid structure; 

1) obtaining the next Ui value and determining, for all existing 

columns, whether the centroid is within the region; 

m) assigning the centroid as the last centroid in the given 

column if the centroid is within the region, or, assigning the centroid 
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as the last centroid in a new column if the centroid is outside the 
region; 

n) repeating steps (l-m) for all centroids; and 
o) calculating an average x-position for each column and 
sorting the columns according to the average x-positions to identify a 
first column, Column 1, a next column, Colmnn 2, and so on to 
Column n. 

35. The algorithm of claim 34 wherein the reference position is an 
upper left comer of the image. 

36. A device readable medium having stored thereon an executable 
instruction in the form of the algorithm of claim 15. 

37. A wavefront measuring apparatus, comprising: 

an illumination component adapted to deliver a small spot of 
light onto a retina; 

an imaging component that can form a centroid image of 
illumination light scattered from the retina and exiting a pupil of an 
eye; 

a detector adapted to acquire the centroid image; 

a processor in operative connection with the detector that can 
execute a centroid displacement calculation to determine wavefront 
measurement infonnatipn; 
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a display component in operative connection with the processor 
that can display a selected wavefront measurement information; and 

means for instracting an online calculation and display of the 
selected wavefront measurement information. 
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